Oscillatory instabilities in d.c. biased quantum dots 
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We consider a 'quantum dot' in the Coulomb blockade regime, subject to an arbitrarily large 
source-drain voltage V. When V is small, quantum dots with odd electron occupation display 
the Kondo effect, giving rise to enhanced conductance. Here we investigate the regime where V is 
increased beyond the Kondo temperature and the Kondo resonance splits into two components. It is 
shown that interference between them results in spontaneous oscillations of the current through the 
dot. The theory predicts the appearance of "Shapiro steps" in the current-voltage characteristics 
of an irradiated quantum dot; these would constitute an experimental signature of the predicted 
effect. 
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Quantum dot devices [jj] display very rich physics in 
and near equilibrium. In particular, experiments Q over 
the past few years have confirmed a decade-old theoret- 
ical prediction that these systems exhibit the Kondo ef- 
fect H : the low temperature formation of a hybridization 
resonance between the dot and the leads to which it is 
connected. This resonance is the result of collective spin 
exchange processes between the leads and the dot, which 
dominate the low temperature physics below a certain 
scale (the 'Kondo temperature'). The Kondo effect is 
manifested as an enhancement of the dot's conductance. 

A bias voltage applied to a quantum dot modifies the 
electron energies in the leads, driving a current through 
the dot. This offers a unique opportunity to study cor- 
related electrons out of equilibrium. Can new kinds of 
non-equilibrium collective phenomena develop in driven 
quantum systems? In the realm of classical physics, there 
are many instances of transitions to new phases in re- 
sponse to a non-equilibrium driving force. An example is 
the Rayleigh-Benard instability, where a temperature dif- 
ference causes the development of convective roll patterns 
that spontaneously select their own spatial and temporal 
frequencies [|J. However, for a driven quantum system 
to develop new collective behavior, it must preserve its 
quantum mechanical coherence. The issue of whether a 
bias voltage preserves the phase coherence of a quantum 
dot is thus a matter of considerable concern. 

In this paper we argue that electron flow through a 
quantum dot modifies, but does not dephase the correla- 
tions of the dot. Central to our arguments is the obser- 
vation that the physics of a quantum dot can be mapped 
onto a one-dimensional problem || , where the electrons 
in each lead are represented by waves traveling in one di- 
rection along one-dimensional conductors, arriving at the 
dot from the left and leaving it to the right (Fig. |l]) . This 
chiral mapping is exact. Electrons scattering off the dot 
do not lose phase information until they suffer subsequent 



inelastic collisions with phonons or magnetic impurities 
in the leads, which, in the chiral mapping picture, lie to 
the 'right' of the dot. Since information travels only from 
left to right, the dephasing of an outgoing electron does 
not affect the spin it has left behind, and thus a current 
of electrons through a quantum dot does not dephase the 
Kondo effect, but acts as a coherent driving force. 
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FIG. 1. (a) A schematic depiction of an incoming electron 
wave A scattering off the quantum dot into the same (B) 
and the opposite (B') leads. Dephasing (indicated by wavy 
crosses) takes place after the scattering event, (b) The map- 
ping of the quantum dot problem to a one dimensional "chi- 
ral" model. Electrons and information travel from left to 
right. Dephasing of electrons at B and B' cannot affect the 
quantum dot. 

This argument is employed to examine the possibility 
of dynamical instabilities in a driven, yet fully coher- 
ent quantum dot. We present calculations which predict 
that, beyond a critical bias voltage V c , the current flow- 
ing through the quantum dot spontaneously acquires an 
oscillatory component whose frequency is a non-linear 
function of the bias voltage. 

The quantum dot system considered here consists es- 
sentially of three parts: the dot itself, and the left and 
right leads, to which it is connected by weak tunneling 
junctions. In the Coulomb blockade regime, the occupa- 
tion of the dot is fixed, and the only allowed processes 
involve flipping the spin of the unpaired electron on the 
dot. In this simple picture, a dot containing an odd num- 
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ber of electrons is a spin-i impurity. As the temperature 
is lowered, spin fluctuations lead to the development of a 
sharp resonance peak in the interacting electron density 
of states: the Abrikosov-Suhl (AS) resonance. This peak 
occurs at an energy tiwi, close to the common Fermi en- 
ergy of the leads, and has a width A ~ Tk, the Kondo 
temperature. (See Fig. ||(a).) 
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FIG. 2. (a) Schematic picture of the equilibrium situation: 
a single Kondo resonance is formed between the electron on 
the dot and electrons in both leads, (b) Quantum dot under a 
DC bias V > V c : Kondo resonances develop at each chemical 
potential, but each resonance still involves the electrons from 
both leads. Interference between the resonances drives the 
oscillatory current. 

In the description which we shall use, the coherence 
associated with the Kondo effect is encoded in terms of 
a time-dependent hybridization field Tp(t) ~ ^ c 0v.adtj 

between the lead electrons and the electron on the dot. 
(The operator cg k(T creates a lead electron with momen- 
tum k and spin a in lead j3 £ {L, R}; d a annihilates a 
dot electron with spin a.) Within the mean field for- 
malism (see below), T^ is replaced by its mean value 
r ee (f>) = Toe^. 

In the presence of a sufficiently large bias voltage V, 
the AS resonance splits into two peaks || , one near each 
Fermi level. (See Fig. |](b).) T(t) therefore becomes a 
sum of two terms, T(t) = Tie iuJlt + T 2 e iuJ2t , where u x and 
L02 are the level positions of the two AS peaks. Quan- 
tum interference between these two terms produces an 
oscillatory current at frequency (cji — 1V2) ~ eV/fi. We 
emphasise that this oscillating current originates as a re- 
sponse to a DC bias: this is to be contrasted with the 
case where the driving voltage is itself time-dependent 

A natural choice for the description of a quantum dot 
coupled to two leads is the Anderson model. We con- 
sider a generalized version of such a model, allowing for 



an asymmetry between transmission and reflection am- 
plitudes (see below). In the Coulomb blockade regime, 
charge fluctuations are virtual, and may be integrated 
out H to obtain an effective spin-exchange Hamiltonian: 

H = H + H R + H T ; (1) 

H = Y1 ( £k ~ ^) C /3ka c /3k<7, 
/3k<r 

Hr = ~ 2~N ^ \_ c La d <y d \c LT + {L^R) , 

Ht = -(1 - ri)— \c Ra d a d\c LT + {L^R) , 

err 

where cl CT = X) c /3k<r- The s P m variables <7,r now take 

k 

integer values from 1 to N, a generalization which we 
discuss below. Here Hq controls the electrons in the 
leads, with /it, — /j_r ee eV, the source-drain voltage. The 
remaining parts of the Hamiltonian describe the spin 
exchange processes which accompany electron reflection 
(H R ) and transmission (Ht) at the dot. 

A crucial ingredient is the asymmetry parameter r\. 
When (Q) is derived directly from the Anderson model, 
one finds that 77 = 0; this reflects a symmetry between 
reflection and transmission processes that will never be 
exactly realized in a physical quantum dot system ||. 
Moreover, it has been proposed |l(| that, at voltages 
above the Kondo temperature, the dot system devel- 
ops a divergent susceptibility to perturbations of the 
form — XO^O, where operator O is defined by 0(t) = 
Y (c Lcr cr aT CR T \ ■ S, where S is the impurity spin. The 

err ^ ' 

inclusion of a small non-zero 77 represents the effect of 
such terms. 

Our analysis of the model ([!]) is based on the 
self-consistent determination of the hybridization fields 
f/3 CT (t) = Cg a (t)d a (t). In an actual physical system, the 
spin index a may take one of two values (f or j), but it 
is useful to consider a general case where a € {1, ... , N}. 
This is because, as N — » 00, the hybridization fields 
r/3 CT (t) behave as well defined semi-classical variables: 

approach deserves particular discussion. By making this 
choice we render the problem exactly solvable, but in 
doing so, we exclude inelastic scattering processes and 
thus emphasize the coherent aspects of the quantum dot 
physics. (Specifically, the lowest order processes that 
produce a lifetime for the rf-fermion are of 0(1/N), and 
therefore disappear in the N — > 00 limit.) The phase co- 
herent physics of the N — 2 Kondo effect is captured by 
the large- N limit pT| . Thus, if an applied voltage acts as 
a phase coherent driving force on the quantum dot, we 
expect the large- TV mean field theory to correctly cap- 
ture any collective behavior that develops at large volt- 
age bias. Support for this approach is provided by re- 
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cent calculations |Tc}], which show that the quantum dot 
system remains a non-perturbative phenomenon at arbi- 
trarily large voltages. However, the ultimate test of the 
procedure must surely lie in experiment. 

Replacing the hybridization fields by their expecta- 
tion values in the large- AT limit results in the mean-field 
Hamiltonian 

H = H Q + J2 [(V£(i)4A + 4cLaV L (t)) + (L - Rj) 

(7 

N\V L +V B \ 2 N\V. 



+ ~ J(2-V) ' JV 
where Vl and Vr are defined by 



V R \ 2 



Vn(t)J 



J JQ.-V) 

J(l-Tj) J 



r«(<) 



(2) 



(3) 



The amplitudes Vl._r(0 must be self-consistently deter- 
mined from the expectation values TL t n(t) at finite volt- 
age V. 

Our analysis uses the Schwinger-Keldysh formalism 
fl2"|| , and from the practical point of view this requires 
an Ansatz for the form of the hybridization fields Vp(t). 
In the present work, two Ansatze are investigated: a two- 
frequency Ansatz ("2F"), 



V = 



Vr 



= Vo(£ H 



(4) 



and a one- frequency Ansatz ("IF"), where the second 
term on the RHS of ([IJ) is absent. The IF solutions cor- 
respond to states with a single AS resonance. The 2F 
solutions correspond to states in which the resonance has 
been split into two components at energies ±w , giving 
rise to an oscillatory current. 

We supplement this analytic approach with numeri- 
cal simulations, using the Heisenberg equations of mo- 
tion to time-evolve the averages ^Cg CT (t)d CT (f)\ and 

( c acr(*) c /3cr(*))j which together describe the state of the 
whole system. The initial state is taken to be a 'high- 
temperature' state, i.e. one with extremely small hy- 
bridization between the dot and the leads. The fields 
Vp(t) are updated at each time step to ensure self- 
consistency. The quantity of interest is the physical cur- 
rent through the dot, J(i) = 2J(1 - rj)Im(Y R (t)Y* L (t)). 

Consider first the IF solutions. These are static solu- 
tions since, with only one AS resonance present, there is 
no possibility of interference: the current through the dot 
is hence time- independent. As long as 77 7^ 0, static solu- 
tions exist for all values of V. As V is increased from zero, 
the level position uj q = does not change, i.e. the applica- 
tion of a small bias neither splits nor shifts the AS peak. 
The width of the peak, A, decreases sharply as one in- 
creases the voltage (Fig. |(b)). At V = V c » 2T K /e, the 
single static solution at frequency lj q bifurcates, and for 



V > V c one finds two degenerate static solutions, with fre- 
quencies ±w , where lj = |^ §(eV/2k B T K ) (Fig. |(a)). 
At large voltages, $>(x) — > 1; thus these degenerate so- 
lutions represent an AS resonance, of width A(?y, V), at- 
tached to the left or right lead. 
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FIG. 3. (a) The unstable frequency u> (V) governing the 
splitting of the AS resonance, in the case n sa 0.37. (b) The 
width A(V) of the AS resonance for the case n » 0.37. 

When V < V c , the numerical simulations support these 
results; however, for V > V c , they indicate that the static 
solutions are unstable to the emergence of a 2F solution 
(||), consisting of components at both frequencies ±uj . 
Quantum interference between these two AS resonances 
contributes an oscillating part to the current I(t) (Fig. ||). 
The oscillations are monochromatic, and their frequency 
is indeed uj a (Fig. ||(a)). 
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FIG. 4. Tunneling current through the dot, I(t), as a func- 
tion of time, for DC voltage V = 3.5 Tk ■ (Here n « 0.46, and 
time is measured in units of the Kondo time, tk = Ti/ksTK-) 

For the oscillating current I(t), one may identify / max 
(the maximum instantaneous current) and J m i n (the min- 
imum); one may thence define the DC and AC parts: 

^DC = 2 (^niin + ^max)i -^AC = 2 C^nax — Imin)- These 

are plotted in Fig. || as functions of V. The appear- 
ance of alternating and direct components to the cur- 
rent is a manifestation of a two-frequency solution, for 
if the hybridization fields have a two frequency form 



T R (t) = rj(-t) = Y+c 



Y*_e ZUJot , then the current 



I{t) = 2J(1 - rf) Im {Y R {t)Y* L {t)) takes the form 
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2 J(1 _ } = M^l + r 2 _] + 2|r+r_ | sin(2a, t + 0), (5) 

where we have written r + r_ = |r + r_|e~ 4 *. Phase co- 
herence between T + and T_ is needed to stabilize the 
phase <fr inside the oscillatory term. 
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FIG. 5. The direct and alternating parts of the current 
through the dot (7dc and Iac) as a function of voltage, V, in 
the case i) » 0.46. 

At small V, the I— V curve is linear, and then peaks 
at eV ~ UbTk- The behavior of the differential conduc- 
tance G(V) = dlnc/dV beyond the peak is novel, as it 
is affected by the sharp transition into the new regime 
where direct and alternating currents occur together. As 
the voltage is increased from zero, G(V) changes sign at 
eV ~ ksTx, exhibits a discontinuity at V = V c , and 
eventually approaches zero from below. 

We conclude with a brief discussion on the detection 
of this proposed oscillatory phenomenon. In close anal- 
ogy to the Josephson effect, we expect that when the 
quantum dot is exposed to radiation, the current-voltage 
profile will develop Shapiro steps [|l3| when the oscilla- 
tory frequency is commensurate with the frequency of 
the incident radiation. We give here a brief derivation of 
the expected effect. 

A modulation in the bias voltage V(t) = V +Vi sin u>i n t 
will lead to a modulation in the separation of the Kondo 
resonances on opposite leads, given approximately by 
2u(t) — 2uj + ^ coauji n t. Since this separation de- 
termines the rate of change of the oscillatory current's 
phase, the current in the irradiated dot will take the form 

/(*) = Ioc + Iac sin ^2 J w{t')dt^j 

= Ibc + Iac sin ^2u> t + sin(w in i) + <p \ . (6) 

This function contains new Fourier components at fre- 
quencies 2oj ± nwj n , where n is an integer, so that when 
2lu is a multiple of LUi n , additional contributions appear 
in the direct current. In a current biased quantum dot, 
we expect this rectification effect to lead to Shapiro steps 



in the current at voltages where u!- m is commensurate 
with the separation of the Kondo resonances 2oj (V). Ex- 
perimentally, these steps should be visible in the range 
Tk < V < U, where U is the charging energy of the 
dot. Observation of such Shapiro steps would constitute 
direct evidence of phase coherent current oscillations in 
the DC biased quantum dot. 
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